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Abstract 

We analyze the role of the nonperturbative input in the Collins, Soper, and 
Sterman (CSS)'s 6-space QCD resummation formalism for Drell-Yan trans- 
verse momentum (Qt) distributions, and investigate the predictive power of 
the CSS formalism. We find that the predictive power of the CSS formal- 
ism has a strong dependence on the collision energy y/S in addition to its 
well-known Q 2 dependence, and the \fS dependence improves the predictive 
power at collider energies. We show that a reliable extrapolation from per- 
turbatively resummed 6-space distributions to the nonperturbative large b 
region is necessary to ensure the correct Qt distributions. By adding power 
corrections to the renormalization group equations in the CSS formalism, we 
derive a new extrapolation formalism. We demonstrate that at collider ener- 
gies, the CSS resummation formalism plus our extrapolation has an excellent 
predictive power for W and Z production at all transverse momenta Qt < Q- 
We also show that the 6-space resummed Qt distributions provide a good 
description of Drell-Yan data at fixed target energies. 

PACS number(s): 12.38.Cy, 12.38.Qk, 14.70.-e 
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I. INTRODUCTION 



The perturbation theory of Quantum Chromodynamics (QCD) has been very successful 
in interpreting and predicting high energy scattering processes. With new data from Fer- 
milab Run II and the LHC in the near future, we expect to test QCD to a new level of 
accuracy, and also expect that a better understanding of QCD will underpin precision tests 
of the Electroweak interactions and particle searches beyond the Standard Model [|l|. As 
pointed out in Ref. [Q], the description of vector and scalar boson production properties, 
in particular their transverse momentum (Qt) distribution, is likely to be one of the most 
intensively investigated topics at both Fermilab and the LHC, especially in the context of 
Higgs searches. It is the purpose of this paper to investigate the predictive power of QCD 
perturbation theory for the transverse momentum distributions of vector and scalar boson 
production in hadronic collisions. 

The production of vector bosons (V = 7*, W ± , and Z) with invariant mass Q at large and 
small Qt in hadronic collisions has been extensively studied |L|]. When Qt ~ Q, effectively, 
there is only one hard momentum scale in the production. Therefore, we expect the fixed- 
order perturbative calculations in power series of ot s to be reliable 0. In this paper, we 
will concentrate on the production of vector bosons at small transverse momenta: Qt < Q, 
where the bulk of the data is. The small Qt region also corresponds to a phase space that 
is most relevant to the hadronic Higgs production. 

When Q T <C Q, the Qt distributions calculated order-by-order in a s in the conventional 
fixed-order perturbation theory receive a large logarithm, hi(Q 2 /Q T ), at every power of a s . 
Even at the leading order in a s , the cross section da/dQ 2 dQ\ contains a term proportional to 
(a s /Q T ) ln(Q 2 /Q T ) coming from the partonic subprocess: q + q — > V(t*, W/Z) +g. Beyond 
the leading order, we can actually get two powers of the logarithm for every power of a s , 
due to soft and collinear gluons emitted by the incoming partons. Therefore, at sufficiently 
small Qt, the convergence of the conventional perturbative expansion in powers of a s is 
impaired, and the logarithms must be resummed. 

Resummation of the large logarithms in QCD can be carried out either in QT-space 
directly, or in the impact parameter, 6-space, which is a Fourier conjugate of the Qr-space. 
It was first shown by Dokshitzer, Diakonov and Troyan (DDT) that in the double leading 
logarithm approximation (DDLA), which corresponds to the phase space where the radiated 
gluons are both soft and collinear with strong ordering in their transverse momenta, the 
dominant contributions in the small Qt region can be resummed into a Sudakov form factor 
0. However, the strong ordering in transverse momenta in the DDLA overly constrains the 
phase space of the emitted soft gluons and ignores the overall momentum conservation. As 
a result, the DDT resummation formalism unphysically suppresses the Qt distributions at 
small Qt 0. By imposing the transverse momentum conservation without assuming the 
strong ordering in transverse momenta of radiated gluons, Parisi and Petronzio introduced 
the 6-space resummation method which allows resummation of some subleading logarithms 
H. By using the renormalization group equation technique, Collins and Soper improved the 
6-space resummation to resum all logarithms as singular as, \n m (Q 2 /Q T )/Q T , as Qt — > 
|| . In the frame work of this renormalization group improved 6-space resummation, Collins, 
Soper, and Sterman (CSS) derived a formalism for the transverse momentum distributions 
of vector boson production in hadronic collisions ||. This formalism, which is often called 
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CSS formalism, can be also applied to hadronic productions of Higgs bosons 0. 

For the Drell-Yan vector boson production in hadronic collisions between hadrons A and 
B, A(P A ) + B(P B ) -> V(Q) +X with V = 7*, W ± , and Z, the CSS resummation formalism 
has the following generic form || 



da(h A + h B ^V + X) 1 



dQ 2 dydQ 2 T (2tt) 



d 2 6 e^' S W(6, Q, x A , x B ) + F (Q T , Q, x A , x B ) , (1) 



where xa = e y ' Q/y/~S and x B = e y Q/y/S with the rapidity y and collision energy v^. In 
Eq. (HI), the term dominates the Qt distributions when Qt <C Q, and the Y term gives 
corrections that are negligible for small Qt, but becomes important when Qt ~ Q- In the 
CSS formalism, the W has the following form || 

Q, x A , x B ) = Wij(b, Q, x A , x B ) <Jij->v(Q) , (2) 

ij 

where a^^viQ) is the lowest order cross section for a pair of quark and antiquark of invariant 
mass Q to annihilate into a vector boson V, and the J2ij r un over all possible quark and 
antiquark flavors that can annihilate into the vector boson at the Born level. In Eq. (|J), 
the Wij(b, Q, xa, %b) is an effective flux to have partons of flavor i and j from the respective 
hadron A and B, and it has the following form |J 

W tJ (b, Q, x A , x B ) = e s ^ W l3 {b, c/b, x A , x B ) , (3) 

where S(b,Q) will be specified later and c is a constant of order one ||@]. The 
Wij(b,c/b,x a,xb) in Eq. (^) depends on only one momentum scale, 1/6, and is perturba- 
tively calculable as long as 1/b is large enough. All large logarithms from ln(l/6 2 ) to \n(Q 2 ) 
in Wij(b,Q,XA,xs) are completely resummed into the exponential factor exp[S(b,Q)]. 

Since the perturbatively resummed Wij(b,Q,XA.,xg) in Eq. (|3]) is only reliable for the 
small b region, an extrapolation to the large b region is necessary in order to complete the 
Fourier transform in Eq. ([I]). In the CSS formalism, a variable 6* and a nonperturbative 
function F NP {b,Q,XA,XB) were introduced 0, 

W css (b, Q, x A , x B ) = W(K, Q, x A , x B ) F NP (b, Q, x A , x B ) , (4) 



where b* — 6/yl + (b/b max ) 2 < b max = 0.5 GeV 1 , and F NP has a Gaussian-like dependence 
on b, F NP ~ exp(— nb 2 ) and the parameter k has some dependence on Q 2 , xa, and xb- 

The predictive power of the CSS 6-space resummation formalism relies on the follow- 
ing criterions: (1) the Fourier transform (or 6-integration) in Eq. ([I]) is dominated by the 
perturbative (or small b) region; and (2) the nonperturbative input F NP has a derived Q 
dependence when Q is large. By using data at some values of Q to fix the nonperturbative 
6-dependence of F NP , the CSS formalism predicts the Qt distributions at different values 
of Q. Existing data are not inconsistent with such a form However, improvements 

are definitely needed for the precision tests of the theory |HJl4] . 



Although the 6-space resummation formalism has been successful in interpreting existing 
data, it was argued []I|,|Ilf that the formalism has many drawbacks associated with working 
in impact parameter space. As listed in Ref. M, the first is the difficulty of matching the 
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resummed and fixed-order predictions; and the second is to know the quantitative difference 
between the prediction and the fitting because of the introduction of a nonperturbative 
F NP . In viewing of these difficulties, major efforts have been devoted to resum the large 
logarithms directly in Q T -space |l|,0 . 

In this paper, we argue and demonstrate that both of these drawbacks can be overcomed 
||. Since there is no preferred transverse direction, the two-dimensional Fourier transform 
in Eq. ([!]) can be reduced into an one-dimensional integration over 6 = |6| weighted by the 
Bessel function J Q {Q T b) We find that by using an integral form for the Bessel function, 
the 6-space resummation formalism works smoothly for all Qt < Q. Because of the smooth 
resummed Qt distributions, the transition (or switch over) to the fixed order perturbative 



calculations at large Qt becomes less ambiguous WA 



In addition, we find that the predictive power of the 6-space resummation formalism 
strongly depends on the collision energy \/~S in addition to its well-known Q 2 dependence. 
Because of the steep evolution of parton distributions at small x, the V~S dependence of 
the W in Eq. (|l|) significantly improves the predictive power of the 6-space resummation 
formalism at collider energies. We demonstrate that the 6-space resummation formalism 
has an excellent predictive power for the Drell-Yan Qt distributions as long as Q 2 is large 
and/or xa and xb are small. 

To quantitatively separate the QCD prediction from the parameter fitting, we derive a 
new functional form in Eq. (^) to extrapolate the perturbatively calculated 6-space distri- 
bution W(b, Q, xa, xb) to the large b region. Our extrapolation is derived by adding power 
corrections to the evolution and renormalization group equations in the CSS resummation 
formalism. Our extrapolation preserves the predictive power of perturbative calculations at 
small b, while it provides clear physical interpretations for all 6-dependence in the large b 
region. We find that the CSS 6-space resummation formalism plus our extrapolation gives 
an excellent description of the data on the Drell-Yan Qt distributions at both collider and 
fixed target energies. 

The rest of this paper is organized as follows. In Sec. [II], we briefly review the CSS 6-space 
resummation formalism for the Drell-Yan transverse momentum distributions. We show that 
the predictive power of the 6-space resummation formalism has a significant v^S dependence. 
In Sec. |TTTj, we quantitatively analyze the role of the nonperturbative input at large b in the 
6-space resummation formalism. By adding power corrections to the evolution equation of 
Wij and power corrections to the renormalization group equations of corresponding evolution 
kernels, we derive a new functional form to extrapolate the perturbatively calculated to 
the large b region. This new functional form clearly separates the perturbative predictions 
in the small b region from the nonperturbative physics at large b. Finally, in Sec. ffVL we 
numerically compare the 6-space resummed Qt distributions with experimental data. Our 
conclusions are also given in Sec. [TV[ 



II. THE COLLINS-SOPER-STERMAN FORMALISM 

In this section, we briefly review the CSS 6-space resummation formalism for the Drell- 
Yan transverse momentum distributions. We show that the predictive power of the 6-space 
resummation formalism has a strong dependence on the collision energy v^S" in addition 
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to its well-known Q 2 dependence. We demonstrate that the y/S dependence significantly 
improves the predictive power of the 6-space resummation formalism at collider energies. 

It was shown that for b <C 1/Aqcd, the W(b, Q, xa, %b) is directly related to the sin- 
gular parts of the Drell-Yan Qt distribution as Qt — > 0. More precisely, the W(b, Q, xa, %b) 
includes all singular terms like S 2 (Qt) and [\n m (Q 2 /Q T )/Q T ] TCg with m > 0. The terms that 
are less singular as Qt — > are included in the Y term in Eq. (JJ). The QCD resummation 
of the large logarithms in the CSS formalism is achieved by solving the evolution equation 
for the W, 



ij 



d ~ 
giri Q 2 w ij( b ,Q, x A,x B ) = [K(bfj,,a s (fj,)) + G(QI >, a s (^))] ^(6, x A , x B ) 

and corresponding renormalization group equations for the kernels K and G, 



d 



<91n/z 

d 



2 K(bfi,a s {n)) = --7 K (a s (//)) 



d In fi 



(5) 

(6) 
(7) 



The anomalous dimension 7x(a s (yu)) = J2 n =i 1 < k\ q 's{^) l^) n 111 Eqs. @ and (0) are pertur- 
batively calculable f|. The renormalization group equations for the K and G in Eqs. (f|) 
and (0) ensure the correct renormalization scale dependence, d/d\iafi 2 [W(b, Q, xa, %b)] = 0. 
The solution given in Eq. (H) corresponds to solving the evolution equation in Eq. (JS)) from 
ln(c 2 /6 2 ) to ln(Q 2 ), and solving the renormalization group equations in Eqs. (]^) and (^) 
from ln(c 2 /6 2 ) to ln(/i 2 ) and from ln(Q 2 ) to ln(/x 2 ), respectively. 

Integrating Eq. (g) over ln(/i 2 ) from ln(c 2 /6 2 ) to ln(/i 2 ), and Eq. (|7|) from ln(Q 2 ) to ln(/i 2 ), 
one derives 



K(bfj,, a s (fj,)) + G(Q/fi, a s (n)) 



Q 2 

c 2 /b 2 



da 2 

4r A(a s (fi)) - B(a s (Q)) 
[i 2 



(8) 



where A is a function of 7#(a; s (//)) and K(c,a s (fi)) while B depends on both K(c, a s (Q)) 
and G(l, a s (Q)). The functions A and B do not have large logarithms and have perturbative 



expansions A = £ n=1 A^\a s /-n) n and 5 = £ n=1 £( n )(a a /7r)' 
coefficients in the perturbative expansions are known MM 



respectively. The first two 



A® = C F 
C F 



2 



N 



67 _ 7T \ _ iu 
18 6 / 9 Rn/ 



10. 



S (2) 



--6 F 
2 



2 r 



7T 



4 

17 2 2 
IF ~ 9^ 



12C(3) 



2 



11 



-7T 



193 
~2A 



+ 3C(3) 



(9) 



where iV = 3 for SU(3) color, C F = (N 2 - 1)/2N = 4/3, T fl = 1/2, and n F is the number 
of active quark flavors. The functions A and B given in Eq. @ are derived from the 
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general expressions in Ref. || with the following choices for the renormalization constants: 
G\ = c — 2e _TB and C 2 = 1, where 7^ « 0.577 is the Euler's constant. 

Substituting Eq. (Q) into Eq. (||), and integrating over ln(Q 2 ) from ln(c 2 /6 2 ) to ln(Q 2 ), 
one obtains the Wij given in Eq. @ with 



5(6, g) 



2 ^ 

c 2 /6 2 /I 2 



In 



'Q 2 



i4(a s (/x)) + £(a s (//)) 



(10) 



In Eq. (H), all large logarithms from ln(c 2 /6 2 ) to ln(Q 2 ) in Wij(b,Q,XA,XB) are completely 
resummed into the exponential factor exp[S(b,Q)], leaving the Wij(b,r,XA,XB) with only 
one momentum scale 1/6. The Wij(b, t,xa,%b) in Eq. @ is then perturbatively calculable 
when the momentum scale 1/6 is large enough, and is given by l6|,[L2jl 



W i:j (b, -,x A ,x B ) 



c c 

fi/A&A, V = 7) fj/B(xB, A* = 7) 



The functions fa/A and fj/B are the modified parton distributions PJT2]|, 

rl xa 
fi/A&A, H) = X/ / ~F Ci/ a { — , fl)) a /A(C, A 1 ) 



(11) 



(12) 



where ^ a runs over all parton flavors. In Eq. fll2|) , (f) a /A(£, A*) is the normal parton dis- 
tribution for finding a parton of flavor a in hadron A, and Cj/ a = S n=0 W/2( Q! s/ 7r ) n are 
perturbatively calculable coefficient functions for finding a parton i from a parton a. The 
first two coefficients of the Cj/ are available 



Ci/l( z i A* 
C$(^ 



c/6) 
c/6) 

c/6) 



0, 

0. 



C§{z,n = c/b) 



y 2 
T fl z (1 - 2: 



1-^) + 



4 5(1 



(13) 



where z and j represent quark or antiquark flavors and g represents a gluon. The coefficient 
functions given in Eq. (|l3|) are derived from the general functional forms in Ref. |J by setting 
the renormalization constants and the factorization scale as, C\ — c, C2 — 1, and fi = c/b. 

The (Tij->v{Q) in Eq. (|J) is the lowest order cross section for a pair of quark and antiquark 
to annihilate into the vector boson (V = 7*, W^, or Z). For V = 7*, we have 



a 



3S 




(14) 



where is the quark fractional charge and N = 3 for SU(3) color. The aij->y(Q) for 
1/ = W ± or Z can be found in Refs. PJT2|. 

In the CSS resummation formalism, the Y term in Eq. (|]) represents a small correction to 
the Qt distribution when Q T <C Q. But, it dominates the (^-distributions when Q T ~ Q. 
The F-term has a perturbative expansion, 

^ = E n =iy {n) M^)/^) n , and the coefficients 

y( n ) have the following factorized form || 
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a b Jx A tZA JiR CR 



x (^f 1 ) R { :Lv{Qt,Q,x a /Ux b /Uh) , (15) 

where J2 a ,b run over all possible parton flavors and /i represents both the factorization and 
renormalization scale. The R^^y * n ^1- (HID are perturbatively calculable and have the 
same normalization as those introduced in Ref. ||. The leading order terms i^L^* are 
available and are given by Eqs. (2.9) to (2.12) in Ref. 0. For and Z production, one 



needs to change the fractional quark charge ef in the * by corresponding weak coupling 
constants 0. 

Since the Y term does not have the large logarithms and perturbatively calculable, 
the predictive power of the CSS formalism relies on our ability to predict the W term in 
Eq. ([I]). Because the leading power perturbative QCD calculations and the normal parton 
distributions in Eq. (^) are only valid for /i > fi ~ 1 — 2 GeV, the perturbatively calculated 
6-space distribution W(b, Q, xa, %b) in Eq. ([I]) is reliable only if the momentum scale 1/b > 
/io- On the other hand, the Fourier transform in Eq. ([]]) requires a 6-space distribution 
W(b, Q, xa, xb) for b 6 [0,oo). Therefore, the predictive power of the 6-space resummation 
formalism is limited by our inability to calculate the nonperturbative 6-space distribution 
at large b []T|,p~4 . 



However, the 6-space resummation formalism has a remarkable feature that the re- 
summed exponential factor exp[5*(6, Q)\ suppresses the 6-integral when b is larger than 1/Q. 
Therefore, when Q ^> yUo, it is possible that the Fourier transform in Eq. (|I]) is dominated 
by a region of b much smaller than l/fio, and the calculated Qt distributions are insensitive 
to the nonperturbative information at large b. In fact, using the saddle point method, it 
was shown that for a large enough Q, the QCD perturbation theory is valid even at 
Qt = 0, and the Fourier transform in Eq. (|T|) is dominated by an impact parameter of order 



JSP 



1 / Aqcd \ 



Aqcd V Q 



(16) 



where A = 16/(49 — 2nf) w 0.41 for quark flavors rif = 5. From Eq. flTB]), the momentum 
scale corresponding to the saddle point, 1/&sp ; can be well within the perturbative region if 
the value of Q is large enough. Therefore, the predictive power of the 6-space resummation 
formalism is directly related to the numerical value of the vector boson's invariant mass Q 
II- 

For and Z production, we have Q ~ My/ or Mz and the corresponding momentum 
scale from Eq. (|16|) , 1/&sp ~ IOAqcd ~ 2 GeV, which is at the borderline of the predictive 
power of perturbative QCD calculations without introducing the power corrections. In the 
rest of this section, we show that the next-to-leading order corrections to the function S(b, Q) 
reduces the numerical value of the b$p- Furthermore, we show that the numerical value for 
the saddle point has a strong dependence on the collision energy a/5, and the \/S dependence 
can either improves or reduces the predictive power of the 6-space resummation formalism. 

Since there is no preferred transverse direction, the W in Eq. ([]]) is a function of b = \b\, 
and the Fourier transform can be written as 
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-i- J d 2 be^ rb W(b,Q,x A ,x B ) 

1 r°° ~ c 

/ dbbJ (Q T b)e s ^ ^ i a ij ^ v (Q)W ij (b,-,x A ,x B ), (17) 



2vr Jo 



where Jo(z) with z = Qrb is the Bessel function. In deriving Eq. (|17D, we used Eqs. @ and 
(§). The 6gp in Eq. (16) was derived by solving 



±m(be s ^A =0 (18) 

with only for the function S(b,Q). Solving Eq. (|i~8| ) for the saddle point relies on the 
assumption that the 6-dependence in Wijib, c/b, xa, %b) is smooth around &sp- 

However, we find from Eq. (|TT|) that the 6-dependence in Wij(b, c/b, x A , %b) is strongly 
connected to the numerical values of x A and x B , and can be very important for determining 
the saddle point if xa and xb are very small or very large M. Taking into account the 
full 6-dependence of W^ib, c/b, xa, xb), the saddle point for the 6-integration in Eq. ([□]) at 
Qt = is determined by solving the following equation, 

| In (&e*<*») + | In ^(b, -x A ,x B )) = 0. (19) 

\« ' / b=b 

If the Wij(b,c/b,XA,xs) has a weak 6-dependence around bo, the second term in Eq. (|19"D 
can be neglected, and the b « fe SP . From Eq. ([□]), the 6-dependence of Wij(b, c/b, xa, xb) 
is directly proportional to the evolution of the modified parton distributions: 

d ( c 

— In ^ °ij->v{Q) Wij(b, -, x A , x b ) 



1 



~7T~ fi/A(xA, A*) Or —^—f j/B (x B ,Lx) 

dinfi dm fx 



(20) 



where ii = c/b. Since the coefficient functions C in Eq. (J13I) do not have 6-dependence 
at these orders, the evolution of the modified parton distributions in Eq. ( E2TX ) is directly 
proportional to the evolution of normal parton distributions, (d/dmfx) 0i/a(£)AO- Because of 
the steep falling feature of the normal parton distributions when £ increases, the convolution 



over £ in Eq. (|12| ) is dominated by the value of £ ~ x A - Therefore, the evolution of the 
modified parton distributions in Eq. (|20"D is directly proportional to the evolution of normal 
parton distributions, (d/dln n) 0j/a(£, A*) at £ ~ x A . From the DGLAP equation, it is known 
that (d/d\nn)<p(x, fx) is positive (or negative) for x < xo ~ 0.1 (or x > xo), and the evolution 
is very steep when x is far away from xq- Therefore, the second term in Eq. ([19]) should be 
very important when xa and Xb are much smaller than the xq. 

Since Eq. ( Pf) has a saddle point solution, the first term in Eq. fll9"|) is a decreasing 
function of 6, and it vanishes at b = &sp- Because of the minus sign in Eq. (|20|) and the 
fact that the number of small x partons increases when the scale fx increases, we expect the 
second term in Eq. (|1|) to be negative when xa and xb are smaller than the typical Xq, and 
to reduce the numerical value of the saddle point. As a demonstration, let Q = 6 GeV and 
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y/S=1.8 TeV. Using CTEQ4M parton distributions and A QC d(% = 5) = 0.202 GeV |5| 



one derives from Eq. fll6|) that &sp ~ 1.2 GeV -1 , and might conclude that the perturbatively 
resummed Qt distribution at the given values of Q and y/S is not reliable. However, as 
shown in Fig. 0(a), the integrand of the 6-integration in Eq. ( [171) has a nice saddle point 
at bo ~ 0.38 GeV -1 , which is within the perturbative region. This is due to the fact that 
%a ~ %b ~ 0.003 are very small. The second term in Eq. fll9|) is negative and it reduces the 
numerical value of the saddle point, which is clearly shown in Fig. 0(b). The solid line and 
dashed line represent the first and the second term in Eq. fll9|), respectively. Although the 
solid line in Fig. [TJ(b) never crosses zero for b < 1 GeV -1 , which is consistent with the fact 
that &sp ~ 1-2 GeV -1 , the dashed line is negative and it cancels the solid line to give a nice 
saddle point at b = b ~ 0.38 GeV -1 . 

Similar to Fig. [I], we plot the integrand of the 6-integration in Eq. ([17]) for Z production 
at Tevatron and the LHC energies in Figs. § and 0, respectively. In plotting both figures, 
we used CTEQ4M parton distributions and the perturbatively calculated functions: A, B, 
and C to the next-to- leading order, which are listed in Eqs. (0) and (|13|). From Eq. (|16|). 
we estimate 6gp ~ 0.4 GeV -1 for Q = Mz- As shown in Fig. 0(b), the solid line vanishes at 
b pa 0.27 GeV -1 , which indicates that the inclusion of the B^\ and reduces the 
numerical value of the saddle point, b$p. The dashed line in Fig. 0(b), which corresponds 
to the second term in Eq. (|T9|), further reduces the numerical value of the saddle point 
to bo pa 0.24 GeV -1 . At the LHC energy, xa and xb are much smaller. We then expect 
the second term in Eq. ([19|) to be more important, which is clearly shown in Fig. 0(b). 
The dashed line in Fig. 0(b) has a much larger absolute value in comparison with that in 
Fig. 0(b). Consequently, the numerical value of the saddle point is further reduced from 
b pa 0.24 GeV -1 at y/S = 1.8 TeV to b pa 0.13 GeV -1 at VS = 14 TeV, where the 
perturbative QCD calculations should be reliable. In addition, the narrow width of the 
6-distribution shown in Fig. 0(a) ensures that the ^-integration is dominated by b ~ b . In 
conclusion, even at Qt = 0, the perturbative QCD based 6-space resummation formalism is 
valid as long as the collision energy v^S" is large enough. 

When Qt > 0, the Bessel function J (z = Qr&) further suppresses the large b region of 
the b- integration. Because the argument of the Bessel function is proportional to Qt, the 
large b region is more suppressed if Qt is larger. That is, the larger Qt is, the better the 
6-space resummation formalism is expected to work. However, it has been known [0] that the 
6-space resummed Qt distribution from Eq. (|l|) becomes unphysical or even negative when 
Qt is large. For example, a matching between the resummed and fixed-order calculations 
has to take place at Qt ~ 50 GeV for W ± production when these two predictions cross over 
n|. We will address this puzzle in Sec. [TV]. 

In the rest of this section, we investigate the predictive power of the 6-space resummation 
formalism for the Drell-Yan production at fixed target energies (VS < 40 GeV). Most data 
at the fixed target energies have Q G (5, 12) GeV and Qt at a few GeV or less. From 
Eq. (0), we find that bsp is of order 1 GeV -1 or larger. Because of the low collision energy, 
the typical values of xa and xb are larger than the xo- Therefore, the second term in the 
Eq. (|i9[) should be positive, which increases the numerical value of the saddle point. As an 
illustration, instead of =1.8 TeV, we replot all quantities in Fig. |l| at v^S* = 27.4 GeV 
in Fig. f|, which is the collision energy for Fermilab experiment E288 |16| |. As expected, 



the dashed line is now positive and the saddle point is no longer in the perturbative region. 
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In conclusion, at fixed target energies, the perturbatively calculated 6-space distribution 
derived from the CSS resummation formalism is not sufficient to predict the Drell-Yan Qt 
distributions at Qt = 0. A nonperturbative extrapolation to the large b region is necessary. 

When Q T > 0, the Bessel function J (z = Qrb) suppresses the large b region of the 
6-integration, and improves the predictive power of the 6-space resummation formalism at 
the fixed target energies. In Fig. |5|, we plot the integrand of the ^-integration in Eq. ( |17D at 
Q T = 1 GeV (solid line) and Q T = 2 GeV (dashed line) for Q = 6 GeV and VS = 27 A GeV. 
As shown in Fig. [5], the saddle point for the 6-integration moves to the smaller b region as 
Qt increases. Since the saddle points of both curves in Fig. |5| are within the perturbative 
region, one might expect the 6-space resummation formalism to provide a good description 
of the Qt distributions at these energies. However, due to the oscillatory nature of the Bessel 
function, the precise value of the ^-integration depends on the detailed cancellations in the 
large b region. Therefore, the predictive power of the 6-space resummation formalism at the 
fixed target energies is still limited by our knowledge of the nonperturbative information at 
large b. More discussions are given in Sec. [TV. 



III. EXTRAPOLATION TO THE LARGE B REGION 

In this section, we quantitatively analyze the role of the nonperturbative input at large 
b in the 6-space resummation formalism. We first briefly review the extrapolation defined 
in Eq. (f|) and its status in comparison with the existing data. Then, by adding possible 
power corrections to the renormalization group equations in Eqs. (§) and fl7|), we derive a 
new functional form for extrapolating the perturbatively resummed W(b, Q, xa, %b) to the 
large b region. This new functional form clearly separates the perturbative prediction at 
small b from the nonperturbative physics in the large b region. 



A. Extrapolation proposed by Collins, Soper, and Sterman 



As discussed in last section, the perturbatively resummed W(b, Q, xa, %b) m Eq. (H) 
is only reliable for the small b region. An extrapolation of the perturbatively calculated 
W(b, Q, xa,xb) to the large b region is necessary in order to complete the Fourier transform 
in Eq. (]!]). In Ref. ||, CSS proposed the following extrapolation, 

H> css (6, Q, x A , x B ) = E *«-v(Q) Q, xa, x B ) F% p (b, Q, x A , x B ) , (21) 

u 

where 6* was defined following Eq. (f|), and the perturbatively calculated Wij(b,Q,XA,XB) 
are given in Eq. ([|). The nonperturbative input distributions F^ p have the following func- 
tional form 



Fl? p (b, Q, x A , x b ) = exp - \n(Q 2 /Ql) g^b) - g i/A {xA, b) - gj/ B (x B , b) 



(22) 



where the ln(Q 2 ) dependence is a derived result. The functions gi(b), gi/A{ x A,b), and 
9j/B(x B ,b) are nonperturbative, and should go to zero as b — > 0. The predictive power 
of the CSS formalism relies on the derived Q 2 dependence and the universality of the F^ p . 
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Since the low energy Drell-Yan data are sensitive to the large b region, in principle, one can 
use the low Q 2 data to fix the parameters of the nonperturbative F^ p and predict the Qt 
distributions of and Z production at high Q 2 . 

Davis, Webber, and Stirling (DWS) introduced the following form for the nonperturba- 
tive distribution F^ p § 

F»£ s (b, Q, x A , x B ) = exp [-b 2 ( 9l + g 2 \n(Q/2Q ))\ , (23) 

where Q = 2 GeV = l/b max , and g\ and g 2 are flavor independent fitting parameters. 
Without the flavor dependence, the extrapolated CSS formalism defined in Eq. ([H]) reduces 
to that in Eq. (|). With g x = 0.15 GeV 2 and g 2 = 0.4 GeV 2 , DWS found [|] that the CSS 
6-space resummation formalism gives a reasonable description of the Drell-Yan data from 
Fermilab experiment E288 at VS = 27 A GeV pf as well as CERN ISR experiment R209 
at y/S — 62 GeV QI7|. 



In order to incorporate possible ln(r) dependence with r = Q 2 /S = x A x B , Ladinsky and 
Yuan (LY) proposed a modified functional form for the F^ p [[□]] 

F? p (b,Q,x A ,x B ) = exp [-6 2 ( 9l + g 2 ln(Q/2Q )) - b 9l g 3 ln(100x A x B )] . (24) 

An extra parameter g 3 was introduced in LY's parameterization of the nonperturbative 
F^ p . Similar to DWS' parameterization, no flavor dependence was introduced into the 
nonperturbative distribution. With gi = 0.1 lt^ GeV 2 , g 2 = 0.58±q;2 GeV 2 , and g 3 = 
— 0.15±g;i GeV" 1 , LY was able to fit the R209 Drell-Yan data as well as CDF data on 
W and Z production from Fermilab. More recently, Landry, Brock, Ladinsky, and Yuan 
(LBLY) performed a much more extensive global fit to the low energy Drell-Yan data as 



well as high energy W and Z data by using both DWS and LY parameterizations ||13|| . In 
order to fit both the low energy Drell-Yan and the collider W and Z data, LBLY found that 
it is necessary to introduce a large overall normalization uncertainty in order to include the 



low energy Drell-Yan data (in particular, E288 data) into the global fit [13]. LBLY also 
emphasized that the collider data on Z production are very useful in determining the value 
of fitting parameter g 2 . They concluded |T3j that both DWS and LY parameterizations 
with updated parameters result in good global fits, but give measurable differences in Qt 
distributions of Z production at Tevatron. 

Based on our discussions in last section, the Qt distributions of Z production at collider 
energies should not be very sensitive to the nonperturbative physics from the large b region. 
Any significant dependence on the fitting parameters for Z production would cast a doubt 
on the predictive power of the 6-space resummation formalism. To understand the fitting 
parameter dependence of the Z production, we introduce the following ratio, 

W css (b,Q,x A ,x B ) 



R w (b,Q,x A ,x B ) 



W(b,Q,x A ,x B ) 

W(b*,Q,x A ,x B ) f np^ 
W(b,Q,x A ,x B ) 



Q,x A ,x B ), (25) 



where W(b,Q,x A ,x B ) is the perturbatively calculated 6-space distribution given in Eq. (0). 
In deriving the second line in Eq. (|25"D , we used the fact that both DWS and LY param- 
eterizations of the F NP are independent of the parton flavors. Using CTEQ4M parton 
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distributions, we plot in Fig. ||(a) the ratio W(b*, Q, xa, xs)/W(b, Q, function 
of b at Q = Mz and v^S* = 1.8 TeV. In Fig. |6|(b), we plot the nonperturbative distribution 
F NP (b,Q,XA,XB) as a function of b with both DWS parameters (dashed line) and LY pa- 
rameters (solid line). In Fig. |](c), we plot the ratio Rw(b,Q,XA,XB) defined in Eq. (f25|), 
which is effectively equal to a product of Fig. ^|(a) and Fig. |6](b). From Fig. [(], we learn 
that the introduction of the 6* significantly changes the perturbatively calculated 6-space 
distribution within the perturbative region; and in the same region, the function F NP can 
deviate from the unity by as much as 50 percents for some fitting parameters. In order 
to preserve the predictive power of the perturbative calculations, it is important to keep 
the W css (b, Q, xa, %b) consistent with the perturbatively calculated W(b, Q, xa, xb) when 
b < bmax (or Rw(b, Q, xa,xb) ~ 1)- However, we find that a significant fitting parameter de- 
pendence (as much as 20 percents) was introduced by the W css (b, Q, xa,xb) to the 6-space 
distribution within the perturbative region. The same conclusion holds if we plot the curves 
at different energies or use other sets of fitting parameters available for the F NP JT2,|IB|. 



B. Extrapolation with dynamical power corrections 

In order to separate the perturbative prediction in the small b region from the nonpertur- 
bative physics at large 6, we derive a new functional form to extrapolate the perturbatively 
calculated W(b, Q,xa, xb) to the large b region. Our goal is to have an extrapolation that 
preserves the predictive power of perturbative calculations in the small b region and extends 
to the large b region with as much correct physics as we can put in. 

Taking advantage of our early conclusion that heavy boson production at collider ener- 
gies should not be very sensitive to the large b region, we can improve the leading power 
perturbative QCD calculations by studying the behavior of power corrections in the region 
of 6-space where b is not too much larger than b max . The power correction in QCD is a very 
rich and difficult subject itself [|TE|-|2~T|] . In order to define the power corrections, we have to 



identify a nonperturbative momentum scale, A, which should be of order Aqcd- For exam- 



ple, the nonperturbative scale can be the target mass [22] or matrix elements of high twist 



operators |T8|-20|. In addition, we have to distinguish two different types of power correc- 
tions: (1) the power corrections directly to the physical observables (such as cross sections or 
structure functions), and (2) the power corrections to the evolution or renormalization group 
equations. The type-one (or direct) power corrections are always proportional to the power 
of (A/Q) with the physically observed momentum scale Q |]TB}-f2l][j . Therefore, the effect of 
this type of power corrections to physical observables can be neglected when (A/Q) — > 0. 
Similarly, the type-two (or indirect) power corrections are proportional to the power of (A///) 
with evolution or renormalization scale /i. It is important to note that physical observables 
are not directly proportional to the evolution or renormalization group equations, instead, 
they depend on the solutions of these equations. Therefore, physical observables carry the 
effect of the type-two power corrections for all fi G [Qo, Q] and the boundary conditions at 
the scale Qo p5 |. Even when Q is much larger than A, physical observables can still carry 



a large effect of the type-two power corrections through the evolution from Qq to Q |24]. In 



this subsection, we concentrate on the type-two power corrections. Because they contribute 
to the evolution or renormalization group equations, these power corrections should have a 
dynamical origin. 
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When b > b max , we solve the evolution equation in Eq. fl5|) from ln(c 2 /6^ a2 .) to ln(Q 2 ) 
in order to separate the leading power QCD calculations at small b from the large b region. 
Because the scale Q of the evolution equation is chosen to be larger than c/b max , we can 
ignore the explicit 1/Q 2 power corrections to this equation. 

However, the kernel K(bfi, oe s (ji)) of the evolution equation has an explicit 6-dependence, 
we need to add power corrections to its renormalization group equation when b > b max . 
Similarly, we need to add power corrections to the renormalization group equation of the 
kernel G(Q/ fi, a s (/x)) when > b max . Since we are only interested in deriving a functional 
form of the 6-dependence due to power corrections, we will not attempt to derive the exact 
coefficients of the power corrections to the renormalization group equations by going through 
detailed analysis on the mixing of leading and high twist operators [^3j . Instead, we introduce 
some fitting parameters for the size of the possible power corrections. For including only 
the leading power corrections, we modify the renormalization group equations in Eqs. (||) 
and (0) as follows, 







; K(bfi,a a (n)) 



din fj, 2 

d 1 1 

^G(Q/fi,as(n)) = -7Ar(a s (/i)) + — j K 



<91n/i 



(26) 
(27) 



where 7^ is treated as an unknown parameter here, though it should in principle de- 
pend on a s (fi) and the characteristic size of high twist operators. In Eqs. ( p6|) and 
(P7|), we parameterize the leading power corrections in such a way that they preserve 
[djd In fi 2 )Wij(b, Q, xa, x b ) = 0. Since we are interested in the region of b not too much 
larger than b max , we neglect higher power corrections in Eqs. (|2B[) and (27). This approxi- 
mation is going to be tested in Sec. [TV|- 

By integrating Eq. ( 
ln(/i 2 ), we have 



over ln(/i 2 ) from ln(c jb ) to ln(/r) and Eq. ( f27|) from ln(Q ) to 



K(bji, a a (n)) + G(Q/fM, a s (fi)) 



+ 



du 2 

-4- A(a s (fl)) - B(a s (Q)) 



2 IV? 

m ax 



dfi 2 



c7&: 

2 /6 2 



K{c,oc s {l))-K{c,a s {-^-)) 



1 



o7Jf(a«(/i)) + r^lK 
2 /r 



(28) 



where A and B are the same as those defined in Eq. In deriving Eq. (^), we neglected the 
power corrections for the momentum scale between c 2 /&^ aa . and Q 2 , which is consistent with 
neglecting the 1/Q 2 term in the evolution equation of the Wij(b, Q, xa, xb)- Substituting 
the K and G in Eq. (|28|) into the evolution equation in Eq. @ and solving the evolution 
equation over the ln(Q 2 ) from ln(c 2 jb 2 ma3 ) to ln(Q 2 ), we obtain the solution for b > b max , 



W? z (b,Q,x A ,x B ) 



max 5 

Q,x A ,x B )F ij (b,Q,x A ,x B ;b 

max J 



(29) 



where Wij is the leading power perturbative solution given in Eq. (|), and 
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T(b,Q,x A ,x B ;b max ) = ^(b,c/b maX} x A} x B ) 
i 3 v ai max > W- (h rib t a td^ 

vy i] \ u max i u max i -^Ai B j 



F {j (b } Q,x A ,x B ;b max ) 

V'ij \ymax j C/ O max , X^4, Xfj J 

"2 f (tfy _ (52 )Q \ + Jn ( b 2 _ b 

\\ u ) \ u maxJ J ~ c 2 \ u 

fK(c,a s ( C -))-K(c,a s (-^)) 

\ " "max 



x exp - In ' V maa; 



c 2 



ma.r 



, (30) 



The nonperturbative function i^7 P — >■ 1 as b — > & maa: . In deriving Eq. fl30D , we approximate 
the /^-dependence of 7^(0^ (/i)) in the small /i region by \^ K {a s {jj)) « 7 (fi 2 )~ a with constant 
parameters 7 and a, and we expect a to be much less than one. This approximation is to 
mimic a summation of a perturbative series in powers of the running coupling constant, 



(a s (jd)) m , with the scale /x extrapolated into the nonperturbative region [25]. Since the 
[K(c,a 8 (r)) — K(c,oc s (t- £ —))] in Eq. ( |3~0"D depends only on b and b max through a s , we can 
combine it with the first term and treat the power a and the coefficient as fitting parameters. 
Since the 6-dependence of the Wij(b, c/b max , xa, xb) depends on the evolution of parton 
distributions, the ln[Wy(6, c/b max , x A , x B ) /Wij(b max , c/b max , x A , x B )} changes sign when x A 
and x B are larger or smaller than the typical xq. Therefore, we can parameterize the ratio 

%(b,c/b x A ,x B ) ^ f ]n ( X_AXB\ ( y _ 2 ^ 1 
W^b 

max ; iy u rn ax ■ XA,X B ) { \ X Q J \ / j 

with parameters #3 and /3. In principle, the parameters #3 and /5 as well as xq can de- 
pend on parton flavors because of the flavor dependence of parton evolutions. Since parton 



distributions are near saturation at a very small momentum scale [26], we expect both g 3 



and (5 to be very small, and therefore, this term can be neglected in comparison with other 
terms in the F^ p . Consequently, we can neglect the flavor dependence of the F^ p . In 
conclusion, without losing the characteristic features of the 6-dependence in Eq. (|30|) , we 
can reparameterize the F^ p as 

F$i(b, Q, x A , x B ; b max ) = exp j - ln(^pi) [ 9l ((& 2 )« - (b 2 max ) a ) + g 2 (b 2 - b 2 max )] \ (32) 

where the explicit \n(Q 2 ti^^/c 2 ) dependence is derived from the evolution equation in 
Eq. (H). In Eq. fl32"D , the b 2 term represents the leading power corrections to the renor- 
malization group equations of the kernels K and G; and the (b 2 ) a term is a consequence 
of extrapolating the leading power part of the kernels K and G to the small momentum 
scale 1/b with all powers of running coupling constants resummed |25|j. The actual size of 



g 2 signals the size of dynamical power corrections. The choices for the parameters: g±, g 2 , 
and a < 1, will be discussed in Sec. fV|. 

We summarize this subsection by writing down our derived extrapolation 

W'M,z A ,z B ) = \%% Q ' x £'°> . , l^"- (33) 

[ W{b max , Q, x A , x B ) Fqz {b, Q, x A , x B ; b max ) b > b max 

with the W(b, Q, x A , %b) given in Eq. (§), and the function Fq£ specified in Eq. fl32|). Since 
the evolution equation in Eq. (H) and the renormalization group equations in Eqs. (El) and 
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(0) do not include any power corrections, the solution of these equations, W(b, Q, xa, xb) 
in Eq. (|33"D , is valid only for b < b c with b c ~ 0.75 GeV -1 , which was estimated by setting 
ln(l/6g) ~ b\. Therefore, the numerical value of b max in Eq. fl3"3| ) should not be larger than 
the b c in order to be consistent with the approximation used to derive the W(b, Q, Xa, %b)- 



C. Corrections from parton's intrinsic transverse momentum 

In Eq. (|3^), all 6-dependence in the Fqz are dynamical in nature from the way we solve 
the evolution and renormalization group equations. We show in this subsection that there 
could be corrections from Q-independent intrinsic 6-dependence to the 6-space distribution 
W^ z {b,Q,x A ,x B ) inEq. ©. 

For an arbitrary function F(b), we introduce 

W i:j (b,Q,x A ,x B ) = F(b) Wij(b, Q, x A , x B ) , (34) 

and find that both Wij(b,Q,XA,x B ) and Wij(b,Q,XA,x B ) can be solutions of the same 
evolution equation in Eq. (^). In principle, the function F(b) can also have dependence on 
parton flavors i and j. If we use the Wij(b,Q,x A ,x B ) instead of Wij(b, Q, x A , x B ) as our 
solution for the QCD resummed 6-space distribution in Eq. (g), the function F{b) should be 
very close to one when b is small (< b max ). Otherwise, its inclusion will not be consistent 
with the leading power QCD calculation because the dominant physics in the small b region 
has been included in the perturbatively calculated Wij(b,Q,x A ,x B ). But, when b is larger 
than b max , neither the evolution equation nor the renormalization group equations have any 
constraints on the functional form of the F(b), as long as it is not a function of Q. 

Physically, however, we do not have an arbitrary function F(b). The inclusion of 
any function F(b) should have a correct physics origin. In the above derivation of the 
Wij(b, Q, xa, x b ), we added the missing physics (power corrections) to the QCD resumma- 
tion formalism when b > b max . However, we did not include the effect of partons' intrinsic 
transverse momentum, which should appear as a part of the Wij(b, c/b, xa, x b ) — the bound- 
ary condition for the evolution equation. When b is small, the factorized formula for the 
Wij(b,c/b,XA,x B ) in Eq. ([□]) should be reliable. But, when b is larger than b max or 1/6 is 
of the order of partons' intrinsic transverse momentum, the perturbative QCD factorized 
formula requires a so-called "Qt smearing" to be consistent with experimental data ||27|| . A 
Gaussian-like smearing function is often used and does a good job in interpreting the data 
| 27|| . In 6-space, we can include the effect due to partons' nonvanishing intrinsic transverse 



momentum by choosing 

F(b) =e W (-g 2 b 2 ) , (35) 

with a constant <j2, which should be of order Aq CD . Let gi = Aq CD , we estimate F(b) > 0.99 
for b < b max , which gives literally no effect to the perturbative regime. However, when 
b 3> b max , the F(b) is expected to have a sizable effect on the low energy Drell-Yan data. 

To include the corrections due to partons' nonvanishing intrinsic transverse momentum, 
we modify our extrapolation in Eq. (|33D as follows 



^M,*a,xb) = {%% Q ' X *' b) , pNP(hn , . " h 7 b r (36) 

{ W{b max , Q, x A , x B ) F$£ (6, Q, x A , x B ; b max ) b > b max 
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where the perturbatively calculated W(b, Q, xa, Xb) is the same as that in Eq. d33|). The 
modified nonperturbative function Fq% in Eq. (|36D is given by 

F$£(b, Q, x A , x B ; b max ) = exp | - In [<?i {(b 2 ) a - (b 2 max ) a ) + 92 (b 2 - b 2 max ) 

- g 2 (b 2 - b 2 max ) | . (37) 

Although the terms with g 2 and g 2 have the same 6-dependence, they have different physics 
origins. The term with g 2 represents the effect of partons' nonvanishing intrinsic transverse 
momentum. The term with g 2 comes from the dynamical power corrections, and has an 
explicit dependence on Q. In deriving Eq. (^), we neglected the intrinsic transverse mo- 
mentum corrections exp(— g 2 b 2 ) to the W in the perturbative region, which is consistent 
with keeping only the leading power QCD calculations in this region. 

Our extrapolation defined in Eq. ( |3~6"D clearly separates the calculable perturbative region 
from the large b nonperturbative region. In addition, all 6-dependence in Eq. fl3"T| ) have their 
own physics origins. The (b 2 ) "-dependence mimics the summation of the perturbatively 
calculable leading power contributions to the kernels K and G to all orders in the running 
coupling constant ct s (/j,) with the scale \i running into the nonperturbative region The 



6 2 -dependence is a direct consequence of dynamical power corrections to the renormalization 
group equations of the kernels K and G. We did not include power corrections to the 
evolution equation because of our choice of the b max . We believe that when Q 2 is much 
larger than c 2 ' jb 2 max , our extrapolation defined in Eq. (^) should give a good description of 
the 6-dependence in the region not too much larger than b max , which is most relevant to the 
heavy boson production. Uncertainties for not including higher power corrections can be 
tested by studying the sensitivities on the parameters g 2 and b max M. 



D. Extrapolation to low Q 2 

When we apply our extrapolation defined in Eqs. fl36|) and (|3?D to the low Q 2 Drell-Yan 
data, we might need a few modifications. As discussed in Sec. H, the Fourier transform 
from 6-space to Qy-space is sensitive to the large b region when Q 2 and \/S are both small. 
Therefore, the Qt distributions at low energies are much more sensitive to the parameters 
g 2 and b max or even higher power corrections. 

For deriving our extrapolation, we systematically dropped the power corrections of mo- 
mentum scale between 1/6^^. to Q 2 . When Q 2 is small, the leading power perturbative QCD 
calculations receive relatively larger 1/Q 2 corrections, in particular, the type- two power cor- 
rections. For example, a 1/Q 2 term in the evolution equation in Eq. (||) results into a 
type-two power correction to the perturbatively calculated W(b, Q, xa, Xb) in Eq- (0)- The 
size of the type-two power corrections resummed from A 2 /Q 2 to A 2 /(l/b 2 nax ) is proportional 
to the numerical value of b 2 max A 2 , where A 2 is the characteristic scale of the power cor- 
rections. For a larger b max , the perturbatively calculated W(b,Q,XA,XB) receives a larger 
power correction. 

On the other hand, the leading power contributions are resummed from lnfT/fr 2 ^) to 
ln(Q 2 ), and are proportional to the numerical value of ln(Q 2 ). When Q 2 ~ Af| at collider 
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energies, the ln(Q 2 ) is much larger than the fr^^A 2 , and the power corrections to the pertur- 
batively calculated W(b, Q, x^, %b) can be neglected. However, at the fixed target energies, 
both Q 2 and the leading power contributions are much smaller. Therefore, the power correc- 
tions to the perturbatively calculated leading power contributions become relatively more 
important at fixed target energies. 

Therefore, in order to reduce the relative size of the power corrections in the perturbative 
(or small b) region, we need to reduce the numerical value of b max . On the other hand, we 
prefer to keep the b max as large as possible to have more contributions from the perturbative 
region in order to have more prediction than parameter fitting. The b max = 0.5 GeV -1 was 
proposed in Ref. ||. For our numerical results in Sec. [TV], we will test the sensitivities on 
the choices of b max . 

The ln(Q 2 b^^/c 2 ) dependence in our Fqz is a direct consequence of dropping the type- 
two power corrections to the evolution equation in Eq. ([5]). Therefore, in order to preserve 
the ln(Q 2 b^^/c 2 ) dependence in our extrapolation, we expect to require a smaller b max for 
a better description of the low energy Drell-Yan data. 



IV. NUMERICAL RESULTS AND CONCLUSIONS 

In this section, we numerically compare the Qt distributions derived from the 6-space re- 
summation formalism with experimental data from W and Z production at collider energies 
to the low energy Drell-Yan processes. 



A. Numerical accuracy 

One of the potential drawbacks of the 6-space resummation formalism is the difficulty 
of matching the resummed and fixed-order predictions to the Qt distributions at large Qt 
IJ. It was generally believed [Q] that the 6-space resummed Qt distribution from Eq. ([!]) 
becomes unphysical or negative when Qt is large. For example, a matching between the 
resummed and fixed-order calculations has to take place at Qt ~ 50 GeV for W production 
when these two predictions cross over |TJ| . On the other hand, as we discussed in Sec. [IT], we 



expect the predictions derived from the 6-space resummation formalism in Eq. (p]) to work 
better when Qt is larger because (1) the b- integral is dominated by the smaller b region and 
(2) the perturbatively calculated Y term is larger than the resummed W term. We find that 
this puzzle was mainly caused by the lack of the numerical accuracy of the Bessel function 
used to perform the Fourier transform in Eq. flT7|). As we show below, the Qt distributions 
derived from the 6-space resummation formalism are smoothly consistent with data for all 
transverse momenta up to Q. 

Due to the oscillatory nature of the Bessel function, high numerical accuracy of the Jq{z) 
with z = Qt b in Eq. (|17]) is necessary for ensuring an accurate cancellation in large z region 
for a reliable 6-integration. Because z is proportional to Qt, the number of oscillations 
of the Bessel function strongly depends on the value of Qt for the same range of b. For 
example, when b e (0, 2) GeV -1 , Jo(Qrb) crosses zero 0, 6, and 63 times for Qt = 1, 10, and 
100 GeV, respectively. It is clear that numerical accuracy of the Bessel function is extremely 
important for the large Qt region. We noticed that most work published in the literature 



17 



used some kind of asymptotic form to approximate the Bessel function when z = Qrb is 
large. We find that the usage of an asymptotic form for the Bessel function is a major source 
of the uncertainties observed for the large Qt region. Instead of using an asymptotic form, 
we use the following integral form for the Bessel function 

J (z) = - I" cos (z sin(0)) dd . (38) 

7T JO 

The great advantage of using an integral form is that we can control the numerical accuracy 
of the Bessel function by improving the accuracy of the integration in Eq. fl38l) . One can test 
the numerical accuracy of the Fourier transform in Eq. flT7| ) by using functions whose Fourier 
transform can be carried out analytically. In view of the nonperturbative 6-dependence in 
large b region, we used two functions: exp(— ab) and exp(— b 2 /a 2 ) to test the numerical 
accuracy of the 6-integration in Eq. fllTD . The 6-integration for these two functions can 
be carried out analytically. Having the analytical solution in QT-space, we can study the 
convergence and the numerical accuracy of the Fourier transform at different Qt by varying 
the parameter a of these two functions. We find that by using the integral form of the Bessel 
function, the numerical integration over b defined in Eq. (|T7D is very accurate for a wide 
range of a and Qt, and the accuracy is only limited by the precision of variables used in a 
computer programming language. 



B. Qt distributions of W and Z production 

In this subsection, we compare the predictions of the 6-space resummation formalism 
with Fermilab data on W and Z production, and quantitatively demonstrate the excellent 
predictive power of the 6-space resummation formalism at collider energies ||. We show 
that the 6-space resummed Qt distributions are very insensitive to the parameters in the 
nonperturbative Fq% . 

Since we are not interested in the detailed fitting to the data in this paper, we did not 
perform any simulation on final-state cuts to improve the theory curves in the following 
plots. In all plots, CTEQ4M parton distributions are used; and for the Y term in Eq. ([[]), 
we use \\jQ 2 + Qt f° r t ne factorization and renormalization scale jx in Eq. (|15|). For W and 
Z production, a fixed range of the rapidity was integrated and a narrow width approximation 
was used for Q 2 integration [p8fl . 

We test the sensitivities on the parameters in the Fqz by first setting g 2 = and #2 = 
(no "power" corrections). We then fix gi and a in Eq. ( |3?| ) by requiring the first and 
second order derivatives of the W to be continuous at b = b max = 0.5 GeV -1 , and plot 
our predictions (solid lines) to the Qt distributions of Z and W production at Tevatron in 
Figs. and ||, respectively. In Fig 0, we plot the da/dQT of e + e~ pairs as a function of Qt 
at y/S = 1.8 TeV. The data are from CDF Collaboration p9| . Theory curves (Z-only) are 
from Eq. (fj) with the W given in Eq. (0). Same as in Ref. f29| , an overall normalization 



1.09 was used. In Fig. R we plot the da/dQT for W production with the same b max and 



<?2- The data for W production are from DO Collaboration [30]. For the theory curves, we 
integrate the rapidity of W particle from -2 to +2 and set the overall normalization to be 
one. From Figs. [7] and [8|, it is clear that the QCD predictions from the 6-space resummation 
formalism are consistent with the data for all Qt < Q. 
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We notice that the theory curves in Figs. |7| and |8] are slightly below the data at large 
Qt- We believe that it is because we have only the leading order contribution to the Y term 
in Eq. ([!]). At large Qt, the Y term dominates. Similar to the fixed-order perturbative 
calculations, the next-to-leading order contribution will enhance the theoretical predictions 



3l| . In addition, the inclusion of virtual photon channel and its interference with Z-channel 
should shift the peak of the theory curves to a slightly larger Qt value and make the theory 
curves closer to the data ]T2"| , p9"| l . 

We now test the theory curves' sensitivities on the parameters of Fq£. Since Q is fixed 
to the mass of the vector boson and the terms with g 2 and g 2 have the same 6-dependence, 
we can simplify the following discussion by rewriting the Fq% in Eq. fl3?D for W and Z 
production as follows [§], 

F»?(b, Q, x A , x B ; b max ) = exp { - g x ((b 2 ) a ~ (&LJ Q ) " 92 (b 2 - b 2 max ) }. (39) 

We let g 2 in Eq. (^) be a fitting parameter for any given value of b max and fix gi and 
a by the derivatives. Although the fitting prefers g 2 ~ 0.8 GeV 2 , the Qt distributions 
are extremely insensitive to the choices of b max and g 2 . The total \ 2 are very stable for 
bmax e (0.25,0.8) GeV" 1 and g 2 G (0, 1.6) GeV 2 . In Figs. and g, we also plot the theory 
curves (dashed lines) with g 2 = 0.8 GeV 2 (the best fitting value). Non- vanishing g 2 gives 
a small improvement to the Qt distributions at small Qt- We then vary the value of a 
in Eq. ( |39"D by requiring only the first order derivative to be continuous at b = b max , and 
find equally good theoretical predictions, except very mild oscillations in the curves at very 
large Qt due to the Fourier transform of a less smoother 6-space distribution. It is well- 
known that when Qt is larger, any small kink in the 6-space distributions transforms into 
oscillations in the Qt distributions due to more rapid oscillations from the Bessel function 
JaiQrb). The observed insensitivity on b max , g 2 , and a is a clear evidence that the 6-space 
resummation formalism is not sensitive to the power corrections at collider energies. That 
is, at collider energies, a direct extrapolation of the leading power contributions to the large 
b region (the (b 2 ) a term) represents the most relevant 6-dependence at large b. 

To further test the sensitivities on the nonperturbative F NP , we introduce a cutoff b c to 
the 6-integration in Eq. ([FT 



w(b c ,Q T ) = ±- f bc dbbJ (Q T b)W QZ (b,Q,x A ,x B ), (40) 

Z7T JO 

and define the ratio, 

R(b c , Qt) = —77 -77— ■ 41 

w{b c = oo, Qt) 

Because of the oscillatory nature of the Bessel function, we expect the ratio R(b c , Qt) to 
oscillate around one and eventually converge to one as b c increases. In Fig. || we plot the 
ratio R(b c , Qt) as a function of b c for Z production at different values of Qt- In plotting 
the theory curves in Fig. |9], we set = 1.8 TeV, b max = 0.5 GeV -1 , and g 2 = 0.8 GeV 2 . 
For Q T = 0, 5, 10, and 20 GeV, as shown in Fig. ^ the R(b c , Qt) quickly reaches one 
at b c ~ 1/Qt as b c increases. Even at Qt = 0, R{b C) QT) is close to one within a few 
percents at b c ~ 2 GeV -1 . It tells us that the physics of the Qt distribution is dominated 
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by the perturbative (or small 6) region. However, due to the oscillatory nature of the Bessel 
function, the 6-integral for the Fourier transform in Eq. (|l|) converges roughly at a common 
value of 6 C ~ 2 GeV -1 , which is larger than b max . That is, the nonperturbative extrapolation 
to the large b region is necessary to ensure the correct cancellations in the large b region. 
The small dependence on the parameters g 2 and b max shows that our extrapolation defined 
in Eq. (^) catches the most of the physics in this region, and higher power corrections are 
not important. 

In addition, we show quantitative ^-dependence of the 6-space resummed Qt distribu- 
tions by introducing the following ratio, 

d (f) fj \ — d<J ^2) I da(g 2b ) 

K^t, g 2 ) - dQ2dydQ2T I dQ2dydQ 2 T > W 

where da j ' dQ 2 dydQ T is defined in Eq. (JJ), and g 2b = 0.8 GeV 2 is the best fit value for g 2 . 
In Fig. [Li], we plot the ratio R a (QT,g 2 ) in Eq. ((42]) as a function of Qt for Z production 
at the Tevatron energy. The solid line corresponds to g 2 = 2 g 2b , while the dashed line 
has the g 2 = \g 2b - For almost all Qt < Q, the deviation of the 6-space resummed Qt 
distributions at the Tevatron energy is less than two percents when we change g 2 from half 
to twice of the best fit g 2 value. The small irregularity for Qt > 50 GeV is due to the 
fact that the calculated cross section has about half of a percent numerical uncertainty in 
this region, which is of the same order as the deviation. We find that the deviation at the 
LHC energy is much less than one percent. That is, the Fourier transform from 6-space 
to Qr-space is completely determined by the perturbatively calculated 6-space distribution 
plus our extrapolation, which is totally fixed by the continuity at b = b max . Therefore, 
we conclude that the 6-space resummed Qt distributions for vector boson productions at 
collider energies have an excellent predictive power. 

C. Drell-Yan Qt distributions at low energies 

As discussed in Sec. Q, the Drell-Yan Q T distributions at fixed target energies are much 
more sensitive to the nonperturbative input at large 6. The predictive power of the 6-space 
resummation formalism may be limited. To explore the predictive power of the 6-space 
resummed Qt distributions for the low energy Drell-Yan process, we compare the low energy 
Drell-Yan data with the predictions from the CSS 6-space resummation formalism plus the 
extrapolation defined in Eq. (|5E|). We use a subset of available data to fix the parameters in 
the nonperturbative Fq% , and then, we compare the predictions calculated by using these 
parameters with the rest of the data. 

When Qt < 4 GeV, the Bessel function does not have any oscillation in the perturbative 
region (6 < b max ). A less smoother connection of the W(b, Q, x^, xb) at 6 = b max does not 
produce any apparent oscillation for the Qt distribution in this region. Therefore, we do 
not have to require the derivative of the W® z (b, Q,xa,xb) in Eq. (|36|) to be continuous at 
6 = b m ax- In this subsection, we treat g±, g 2 , g 2 , and a in the Fq% as fitting parameters. In 
order to maximize the predictive power of the 6-space resummed Qt distributions, we try 
to fit the data with the least number of parameters, and then add more parameters to see 
the improvements. By going through this fitting procedure, we can learn the importance of 
each parameter and the corresponding physics. 
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We choose the following Fermilab data to fit the parameters in the Fq£: Q G (5, 6) GeV 
and Q G (7,8) GeV from E288 @; and Q G (7,8) GeV and Q G (10.5,11.5) GeV from 



E605 ||32|| . We have a cut on the transverse momentum: Qt < 1.4 GeV, which is the same 



as that used in Ref. |13| . Since the g 2 dependence is very small at collider energies, we try to 
fit these Drell-Yan data with two parameters: g\ and a plus two overall normalizations for 
two experiments. Using the Fqz i n Eq. ( p7|) and b max = 0.5 GeV -1 , we obtain a convergent 
fit with a total x 2 — 78 for 28 data points. The corresponding numerical values of the fitting 
parameters: g\, a, and two overall normalization constants for these two experiments are 
given in Table |. The large \ 2 clearly indicates that power corrections (or terms proportional 
to b 2 in the Fq%) are very important for understanding the low energy data. 

We now include the b 2 dependence of the Fq£ into the fit. In order to use a minimum 
number of fitting parameters, we first combine both g 2 and g 2 terms and define 

g 2 ln(^^j+g 2 = g. (43) 

We approximate the g as a constant fitting parameter. With three fitting parameters: g±, a, 
and g plus two overall normalizations for two experiments, we fit the same 28 data points for 
b m ax = 0.5 GeV -1 . We obtain a much better fit with a total x 2 = 32, and the corresponding 
fitting parameters are given in Table |. Using b max = 0.5 GeV -1 , and the numerical values 
of the fitting parameters in the Table |, we find that the b 2 term in the Fq% is less than 
the (b 2 ) a term for b as large as b ~ 2.4 GeV -1 (~ 2.0 GeV -1 for b max = 0.3 GeV -1 ). 
It indicates that although the power corrections (the b 2 dependence) are important, the 
extrapolation from the leading power contributions (the (b 2 ) a dependence in the Fq%) is 
crucial for understanding the low energy data. 

As we discussed in Sec. [Ill D| , we can test the size of power corrections to the perturbative 
region (6 < b max ) by studying the b max dependence. We find that the fitting parameter 
a is extremely stable when we change the b max , and it prefers a numerical value around 
a ~ 0.15. This can be understood as follows. The parameter a was introduced in Sec. [111B 
to approximate the /^-dependence when we extrapolate the leading power contributions to 



the S(b,Q) into the small fi 2 region JHJ. Therefore, the parameter a should only depend 
on the effective anomalous dimensions (resummed to all orders in the running coupling 
constant, a s (/i)), and should not be sensitive to the numerical value of the b max . We then 
fix a = 0.15 and re- fit those 28 data points with two parameters: g\ and g plus two overall 
normalization constants. 

We find a smooth reduction of the total x 2 when we decrease the b max to as low as 
0.3 GeV -1 . For b max = 0.5, 0.4, and 0.3, we obtain convergent fits with the total x 2 — 32, 
27, and 24, respectively. The corresponding fitting parameters are listed in Table |. Although 
the total x 2 is still very stable when we further reduce the b max , we stop at b max = 0.3 GeV -1 
because it is difficult to distinguish the perturbative prediction from the parameter fitting 
when b max is too small. We confirm from the fitting results listed in Table | that the 
(b 2 ) a term in our extrapolation is very important and it dominates the transition region 
between the perturbative calculation and the nonperturbative extrapolation. We learn that 
the overall normalizations, which are needed to fit the different experimental data sets, have 
a strong dependence on the b max . As shown in Table Q, the overall normalizations are driven 
to the unity as the b max decreases. That is, if we reduce the relative size of the power 
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corrections in the perturbative region {b < b max ) by a reduction of b max , both data sets used 
in the fit are in an excellent agreement with each other. 

In order to show the quantitative size of the power corrections (the b 2 term in the Fq^) 
at the low energy, we replot Fig. [L0| at Q = 6 GeV and y/S = 27. 4 GeV in Fig. [□]. In 
plotting Fig. |TT], we use the ratio R a (QT,92) defined in Eq. (|42| ) with the g 2 dependence 
replaced by the g dependence defined in Eq. (f43|) . The g 2b in Eq. (^) is replaced by the 
best fit value g>&. We choose b max = 0.3 GeV -1 and the corresponding fitting parameters: 
a, gi, and g given in Table 0. The best fit value g& = 0.28 GeV 2 . In Fig. [□], the solid line 
corresponds to g = 2gt while the dashed line has g = \ gb- For Q T from to 2 GeV, the 
6 2 -dependence can change the 6-space resummed Qt distributions by as much as 80 percent 
when the parameter g changes from half to twice of the best fit value. In contrast to the 
two percents deviation at collider energies, as shown in Fig. [H], the large variation shown 
in Fig. |Tl| further confirms that the power corrections are very important for describing the 
Drell-Yan data at the fixed target energies |53 . 

To explore the predictive power of the 6-space resummed Qt distributions at fixed target 
energies, we compare the resummed Qt distributions with existing Drell-Yan data in Figs. [T3, 
|13| ) and |14[ All theory curves in these figures are calculated with our extrapolation defined 
in Eq. ( PUD and b max = 0.3 GeV^ 1 . The nonperturbative extrapolation to the large b region: 
Fqz is defined in Eq. (|37j) with the g 2 and g 2 terms combined as in Eq. ([4^). We fix the 
parameter a to be 0.15, and use the fitted values g\ = 0.92 GeV a and g = 0.28 GeV 2 from 
Table |. Although we only used 28 data points with Qt < 1.4 GeV in our fits for determining 
the two parameters: g\ and g, we plot both theory curves and data for a much enlarged 
phase space in Figs. [12], [1^, and [14]. We plot data with Qt as large as 2 GeV at different 
values of Q in order to explore the predictive power of the theoretical calculations. In Fig. [12|, 
we compare the theoretical calculations with Fermilab E288 data at \/S = 27 A GeV |L6 



The theory curves are multiplied by an overall normalization constant Ne288 — 0.97 as listed 
in Table |, which is different from what was found in Ref. (T^|. From top to bottom in 
Fig. [TJ, the four curves along with four data sets correspond to Q G (5, 6) GeV, (6, 7) GeV, 
(7, 8) GeV, and (8, 9) GeV, respectively. Two of the four data sets: Q G (5, 6) GeV and 
(7, 8) GeV with Qt < 1.4 GeV were used in our fitting. That is, only 14 of the 40 data 
points plotted in Fig. [12] were used in the fitting. Clearly, for all Qt up to 2 GeV, the 6-space 
resummed Qt distributions are in excellent agreement with the data from E288. In Fig. [13|, 
we plot the resummed Qt distributions along with Fermilab E605 data fl3^[ . The overall 
normalization constant for E605 is one. The exactly same values of g\ and g are used for 
calculating the theory curves in Figs. and [13]. In Fig. from top to bottom, the four 
curves along with four data sets correspond to Q G (7, 8) GeV, (8, 9) GeV, (10.5, 11.5) GeV, 
and (13.5, 18.0) GeV, respectively. Similar to E288 data, only two of the four data sets: 
Q G (7, 8) GeV and (10.5, 11.5) GeV with Q T < 1.4 GeV were used in our fitting. Although 
only 14 of the 40 data points in Fig. were used in the fitting, the 6-space resummed 
Qt distributions are in a good agreement with all 40 points, except a few points with 
Q G (7, 8) GeV. Actually, seven of the ten data points in this set with Q G (7, 8) GeV were 
used in our fitting. However, because of the relatively large error bars, these data points did 
not have enough weight in the fitting. Nevertheless, the theory curves calculated with two 
fitting parameters: g\ and g give a very good description of the low energy Drell-Yan data 
in Figs. ^ and [T3|. In particular, the overall normalization constants for both experiments 
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13], the 5-space resummed Qt distributions have to multiply 



are very close to the unity. 

As pointed out in Ref. 
a large overall normalization constant in order to be consistent with the Fermilab E772 
data. In Fig. [14|, we plot the 6-space resummed Qt distributions along with the E772 data 
|34| . Three data sets from top to bottom correspond to Q G (5,6) GeV, (8,9) GeV, and 
(11, 12) GeV, respectively. The theory curves are calculated with the same formula and 
parameters used to calculate the curves in Figs. |T2| and [H| except a much larger overall 
normalization constant Ne 772 = 1-6, which is consistent with what was found in Ref. 
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Although none of these data points in Fig. [TJ were used in our fitting, the theory curves 
describe the data well. Without plotting another figure, we state that the 6-space resummed 



Qt distributions are also consistent with the R209 data |T7| . 

In order to get a feeling on the relative size of the contributions from the partons' intrinsic 
transverse momentum and the dynamical power corrections, we separate the #2 and <j2 terms 
in Eq. (^). With b max = 0.3 GeV -1 and a fixed a = 0.15, we perform a new fit with three 
fitting parameters: g\, g-i, and §2- For the same 28 data points, we did not get much 
improvement in the total \ 2 ■ This result should be expected because (1) we already have 
an excellent fit with a x 2 /d-o.f ~ 1 when gi and g~2 terms are combined, and (2) the range 
of Q values in our 28 data points are limited. Nevertheless, the fitting result indicates that 
the g2 term and the #2 term are roughly equal. That is, the effect of partons' intrinsic 
transverse momentum, which is Q-independent, i s as important as that of the dynamical 
power corrections at the fixed target energies. 

From all three figures in Figs. |l|, [13|, and [L4], we conclude that the CSS resummation 



formalism in combination with our derived extrapolation provides a very good description 
of the Drell-Yan Qt distributions at the fixed target energies. Using only 28 data points 
with Qt < 1.4 GeV to fix the parameters: g\ and g, our numerical results are consistent 
with over 100 data points from three experiments. In particular, except the E772 data, the 
overall normalization constants between the theory and the data are extremely close to the 
unity. 

To conclude this subsection, we demonstrate the convergence and stability of the Fourier 
transform at the fixed target energies by plotting the ratio R(b c , Qt) as a function of b c 
in Fig. [L5|. The ratio R(b c ,Q T ), defined in Eq. (0), are evaluated at Q = 6 GeV and 
\fS = 27 A GeV. For calculating the resummed 6-space distribution W QZ {b 1 Q, xa-, %b), we 
choose b max = 0.3 GeV -1 and the corresponding fitting parameters listed in Table |, which 



are the same as those used for plotting theory curves in Figs. O, O, and HI. In Fig. [TR 
we plot the ratio R(b c , Qt) for four different Qt values: 0, 0.5, 1.0, and 2.0 GeV. The ratio 
R(P c ,Qt) quickly converges to one at b c ~ 3 GeV -1 . By comparing Fig. [] and Fig. |TR we 
conclude that the Fourier transform for the 6-space resummed Qt distributions converges 
at b c ~ a few GeV -1 , and therefore, higher power corrections are not very important. On 
the other hand, because the Fourier transform converges in the nonperturbative region, a 
reliable extrapolation from the perturbatively calculated 6-space distribution at small b to 
the nonperturbative large b region is necessary and crucial in order to be consistent with 
experimental data. From the consistency shown in Figs. [12|, 13[ and |14, we further conclude 
that the functional dependence on the impact parameter b, introduced in our Fq^ , catches 
the correct physics. 
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D. Conclusions 



In conclusion, we have quantitatively investigated the role of the nonperturbative input 
in the CSS 6-space QCD resummation formalism for the Drell-Yan Qt distributions at both 
collider and fixed target energies. We find that the predictive power of the CSS resummation 
formalism has a strong dependence on the collision energy y~S. The y/S dependence is a 
consequence of the steep evolution of parton distributions at small x, and it significantly 
improves the predictive power of the CSS formalism at collider energies, in particular, at the 
LHC energy. We show that although the resummed Qt distributions are mostly determined 
by the perturbatively calculated 6-space distributions at small b, a reliable extrapolation 
to the nonperturbative large b region is necessary to ensure the correct cancellations of 
the 6-integration when b > 1/Qt- By adding power corrections to the renormalization 
group equations in the CSS resummation formalism, we derive a new functional form in 
Eq. ( |5BD to extrapolate the perturbatively resummed 6-space distributions to the large b 
region. We demonstrate that at collider energies, the CSS resummation formalism with our 
extrapolation has an excellent predictive power for Qt distributions of W and Z production 
for Qt as large as Q. Because of the smooth resummed Qt distributions, the matching 
between the resummed and fixed-order calculations at large Qt is less ambiguous. 

In this paper, we explicitly show that the power corrections are very important for the 
resummed Qt distributions at the fixed target energies. With only two parameters: gi 
and g, obtained by fitting 28 data points, the calculated Qt distributions using the CSS 
formalism plus our extrapolation are in a good agreement with all existing data. In our 
derived extrapolation in Eq. (|6|), the resummed 6-space distributions for b < b max do not 
include perturbative power corrections, while the nonperturbative extrapolation have both 
leading power contributions as well as the power corrections. Therefore, by choosing a 
smaller b max , we effectively move the power corrections to the relatively small b region. We 
find that by reducing b max , the overall normalization constants for E288 and E605 data sets 
are driven to the unity. This result not only shows a good consistency between different 
experiments, but also tells us that the power corrections are very important for describing 
the low energy Drell-Yan transverse momentum distributions. 

Finally, we argue that the CSS 6-space resummation formalism should provide a reliable 
prediction for the Higgs production at the LHC energy. At the Vo = 14 TeV, we expect the 
partonic subprocess: g + g ^ H + X to dominate the Higgs production when Higgs mass 
Mh ~ 115 GeV. From the fact that xa ~ %b ~ 0.008 are small, and gluon distribution has 
a steeper evolution than quark distribution at small x, we expect the VS dependence of the 
W(b, Q, xa, xb) to move the saddle point bo to a value much smaller than the bsp = 0.37 
estimated by using Eq. (jT6|), and most likely, even smaller than 0.13 shown in Fig. ^(b) for Z 
production at the LHC energy. Therefore, the 6-space QCD resummation formalism should 
be valid for predicting the Qt distribution of hadronic production of Higgs bosons at the 
LHC energy. 
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FIG. 1. (a) Integrand of the 6-integration in Eq. ( |17| ) at Qt = and Q = 6 GeV as a function 
of b with an arbitrary normalization at Tevatron energy; (b) the first (solid) and second (dashed) 
terms in Eq. (FL9I) as a function of b at the same Q and y/S. 
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FIG. 2. (a) Integrand of the 6-integration in Eq. ( |17| ) at Qt = and Q = Mz as a function 
of b with an arbitrary normalization at Tevatron energy (y/~S = 1.8 TeV); (b) the first (solid) and 
second (dashed) terms in Eq. ( |T^ ) as a function of b at the same Q and \/S. 
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FIG. 3. (a) Integrand of the 6-integration in Eq. ( |17| ) at Qt = and Q = Mz as a function 
of b with an arbitrary normalization at the LHC energy (y/S = U TeV); (b) the first (solid) and 
second (dashed) terms in Eq. ( |T^ ) as a function of b at the same Q and V^- 
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FIG. 4. (a) Integrand of the 6-integration in Eq. (|l7| ) at Qt = and Q = 6 GeV as a function 
of b with an arbitrary normalization at E288 energy (y/S = 27 A GeV); (b) the first (solid) and 
second (dashed) terms in Eq. ( |T^ ) as a function of b at the same Q and y/S. 
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FIG. 5. Integrand of the 6-integration in Eq. (|l7| ) at Qt = 1 GeV (solid line) and Qt = 2 GeV 
(dashed line) as a function of b with an arbitrary normalization. The Q and are the same as 
those in Fig. |(a). 
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FIG. 6. (a) Ratio W(b*, Q, xa, xb)/W(J), Q, xa, %b) as a function of b; (b) F NP (b,Q,XA,XB) 
as a function of b; (c) Ratio Rw defined in Eq. ( |25| ) as a function of b. All plots have Q = Mz and 
VS = 1.8 TeV. 
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FIG. 7. Comparison between the 6-space resummed Qt distribution and CDF data [29|. The 
inset shows the Qt < 20 GeV region. 
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FIG. 8. Comparison between the 6-space resummed Qt distribution and DO data |^0|. The 
inset shows the Qt < 20 GeV region. 
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FIG. 9. Ratio R(b c ,Q T ) in Eq. @ function of b c for Z production at Tevatron energy: 
Q T = GeV (solid), Q T = 5 GeV (dashed) in the top plot; Q T = 10 GeV (solid), and Q T = 20 GeV 
(dashed) in the bottom plot. 
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FIG. 10. Ratio R ct {Qt,92) in Eq. (g) at Q = M z and y/S = 1.8 TeV as a function of Q T : 
92 = 292 b (solid) and g 2 = (l/2)g 2b (dashed). 
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FIG. 11. Same as that in Fig. 10 at Q = 6 GeV and E288 collision energy. 
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FIG. 12. Comparison between the 6-space resummed Qt distribution and Fermilab E288 data 
[16]. The overall normalization for the theory curves: Ne288 = 0.97. 
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FIG. 13. Comparison between the 6-space resummed Qt distribution and Fermilab E605 data 



[32]. The overall normalization for the theory curves: Ne605 = 1-0. 
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FIG. 14. Comparison between the 6-space resummed Qt distribution and Fermilab E772 data 
[34]. The overall normalization for the theory curves: NE772 = 1-6. 
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FIG. 15. Ratio R(b c , Qt) as a function of b c for Drell-Yan production at Q = 6 GeV and 
VS = 27 A GeV: Q T = GeV (solid), Q T = 0.5 GeV (dashed) in the top plot; Q T = 1 GeV 
(solid), and Qt = 2 GeV (dashed) in the bottom plot. 
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TABLES 



bmax 


x 2 


a 


9i 


9 


Ne288 


Ne605 


(GeV- 1 ) 






(GeV Q ) 


(GeV 2 ) 






0.5 


78 


0.65 


0.4 





0.85 


0.9 


0.5 


32 


0.15 


1.14 


0.19 


0.88 


0.93 


0.4 


27 


0.15 


1.06 


0.23 


0.93 


0.98 


0.3 


24 


0.15 


0.92 


0.28 


0.97 


1.00 



TABLE I. Nonperturbative parameters in the Fqz obtained by fitting 28 data points on 
Drell-Yan Qt distributions at the fixed target energies. 
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